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Abstract 

The electron correlations in narrow energy bands are examined in 
framework of the Hubbard model. The single-particle Green func- 
tion and energy spectrum are obtained in paramagnetic state at half- 
filling by means of new two-pole approximation. In the ground state 
analitical expressions for the energy gap, polar states concentration 
and energy of the system are found. Metal-insulator transitions in 
the model at change of bandwidth or temperature are investigated. 
The obtained results are used for interpretation of some experimental 
data in narrow-band materials. PACS 71.28.+d, 71.27.+a, 71.10.Fd, 
71.30.+h 

1 Introduction 

Among the observed in narrow-bands materials metal-insulator transi- 
tions (MIT) the significant interest is attracted to the transitions from para- 
magnetic metal state to paramagnetic insulator state at increase of temper- 
ature which exhibit the systems NiS^^Se^ (Vi__Cr_) 2 03 [Q, £§] and 
Yx_a.Caa.TiO3 0, in these systems the paramagnetic insulator - paramag- 
netic metal transitions under external pressure are observed also. There are 
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reasons to believe that noted transitions are the consequences of electron- 
electron interactions and can be described within the framework of Hubbard 
model 0| . 

Hubbard model is the simplest model describing MIT in materials with 
narrow energy bands. This model describes a single non-degenerate band of 
electrons with the local Coulomb interaction. The model Hamiltonian con- 
tains two energy parameters: the hopping integral of an electron from one 
site to another and the intraatomic Coulomb repulsion of two electrons of 
the opposite spins. This model is used intensively (for recent reviews see 
Refs. |T|, p[-[P!D[) in order to describe the peculiarities of physical properties 
of narrow-band materials; in this connection two-pole approaches are at- 
tractive. The two-pole approaches in the Hubbard model and the Hubbard 
bands conception (being the consequence of two-pole approximation) have 
been useful for understanding of the peculiarities of electric and magnetic 
properties of narrow-band materials [pi |10[. However within the framework 
of two-pole approaches there are series of issues, in particular the problem of 
metal-insulator transition description |10|, |H|-f2T|j. 

is 



In the present paper recently proposed two-pole approximation [11 



used to study effects of electron correlations in the Hubbard model. The 
single particle Green function and energy spectrum are obtained. In the 
ground state analitical dependences of energy gap, polar states (doublons 
or holes) concentration and energy of the system on model parameters are 
found. Dependences for energy gap, polar states concentration on tempera- 
ture are calculated. The obtained results are compared with corresponding 
results of other approximations and are used for the interpretation of some 
experimental data. In particular, the observable transitions from an insu- 
lating state to a metallic one at increase of bandwidth, and from a metallic 
state to an insulating one at increasing temperature are explained. 



2 Single-particle Green function and energy 
spectrum 

The Hubbard Hamiltonian [1] in terms of transition-operators of z'-site from 
state \l) to state \k) Xf 1 is written as 

H = H + H 1 + H[; (1) 
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H = -^(x-+xf) + uJ2xf, 

icr i 

h x = J2 taixfx^ + xfxf), 

H[= Uj (X?°XJ 2 - XfXf + h.c. 



(2) 
(3) 
(4) 



where /x is the chemical potential, U is the intra-atomic Coulomb repulsion, 
tij is the nearest-neighbor hopping integral, X* is the operator of number 
of |fc)-states on z'-site; a denotes spin of an electron {a =|, t) an d 0" denotes 
the projection of electron spin opposite to a; H describes system in the 
atomic limit, Hi describes electron hoppings between single occupied sites 
and empty sites (holes) (the first sum in Hi - processes forming "h-band") 
and electron hoppings between doubly occupied sites (doublons) and single 
occupied sites (the second sum in Hi - processes forming "d-band"). H[ 
describes "hybridization" between the "h-band" and "d-band" (the processes 
of pair creation and annihilation of holes and doublons). 

The single-particle Green function is written in operators as 

((«ptI4» = «*JW» - - + «* P 0T I*J°». (5) 

The functions ((X^ 2 \X 21 )) and ((X^\X 21 )) satisfy the equations 



{E + n-U)((X^\X^)) 



(£ + A*)«*JW» 



^(Xl+Xl) + (([x^Hi]_\X^)) 
(([xi 2 ,H[]_\X^)), (6) 

((\xy,Hi]jx?)) + {{\xy t Hi]jx>,% 



with [A, B]_ = AB — BA. To obtain the closed system of equations we apply 
new two-pole approximation, proposed in work ]TT[. Suppose in Eq. ([]) that 



0] 

j ' 



Xj; 2 ,Hi]_=J2Kpj)Xj 2 , 



(7) 



where e(pj) and e(pj) are non-operator expressions which we calculate using 
the method of work At electron concentration n—1 in a paramagnetic 
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state we have 



e(pj) = (1 - 2d)t pj , 
e(pj) = (l-2d)t pj , 



(8) 



with d = (Xp) being the concentration of doublons. 

Let us take into account the functions ((Xj ) 2 ,H[ \X^)) and 

(( X®\ H[ \Xg^)) in the mean-field approximation: 



i,i^p 

({xxx?\x?))\ ~ - E M^i + Aj>«jcn^» 

i,i^p 

- E + *, T )*W>> + 

i,i^p 

<<*pWW>>] --E M* P ° + ^><(^W>>; 



in this way we neglect the processes describing the "inter-band" hoppings 
of electrons which are connected with spin turning over and "inter-band" 
hoppings with creation or annihilation of two electrons on the same site. 
So we obtain the closed system of equations 

{E-IJL + U)((X?\X?)) - E e(pi)«Xf \X?)) + (Xi + Xl) £ M<xf 



2tt 



<W 10 ) 



(E - vl){{X*\X?)) - E^)((^ 0T l* 2i >> + <* P ° + * P T > E = o. 

After the Fourier transformation we obtain solutions of system of Eqs. ( P~0| ) : 



2tt 



+ 



S-£7 h (k) E-E d (k) 



(11) 
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1 / _ U-e(k) + i(k) \ Rl _, A 
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x f 1 1 ^ (12 ) 

£ d (k) - E h (k) \E - E h (k) E - E d (k) J ' 1 ] 
Here £(k) is the hopping intergral in k— representation and 

Eh{k) = _,4 + i«±iw 

1 



[U - e(k) + e(k)]» + (X° + X p T >(Xi + AJ)(t(k))«, (13) 



^ ^ f 7 e(k) + e(k) 

1 



[i7-e(k) + e(k)P+(X0 + Xj)(Xi + X2)(£(k))2 (14) 

are the energies of electron in lower ("hole") and upper ("doublon") sub- 
bands, respectively; e(k) and e(k) are the Fourier components of e(pj) and 

Kpj)- 

Analogous procedure gives for functions ((X^ 2 \Xj )) and ((X^|Xj )) the 
following expressions: 

((xf\x!°)) t = ((x'»\xv)) k , 

Finally, in k-representation single-particle Green function (|5|) we obtain 

^ = ^ ( E _^ fc(k) + E -E d (k)) ' (16) 

, If (Ci ~ C 2 )(t/ - e(k) + e(k)) + 4t(k)dC 2 \ 
k 2^ E d (k)-E h (k) 
B k = l- A k , 



where C x = (X° p + X}), C 2 = (X 2 p + 
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In the important for an investigation of metal-insulator transition case 
= 1 in a paramagnetic state ((Xj) = (X};)) single-particle Green func- 



tion ([IB]) has the form 




-If A - + B - ) (17) 
f(k) 



+ (t(k)) 

1-A k , 

where single-particle energy spectrum is 
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^(k) = (l-2d)t(k)--^ 2 + (t(k)) 2 , 



£,(k) = (1 - 2d)t(k) + - ^ + (t(k)) 2 (18) 



(here we took into account that ji = ~ for n = 1). 



Single-particle Green function fll7|) and energy spectrum are exact in 



the band and atomic limits. It is worthwile to note, that in distinction from 
the results of two-pole approximations of Hubbard and Ikeda, Larsen, 
Mattuck [TJ[] the energy spectrum (|Tj|) depends on polar states concentra- 
tion (thus on temperature). In distinction from approximations based on 



ideology of Roth |T5j (in this connection see also Refs. [IS|-||21||) the energy 
spectrum fll8f ) describes metal-insulator transition. Energy spectrum which 
describes metal- insulator transition was earlier obtained in work fOl. Ex- 



pressions ( |I8D differs from the respective expressions in work ]T3[] by presence 



of term y U 2 + t 2 (k) instead of y £/ 2 + 4<i 2 t 2 (k). This leads to the series of 
distinctions between results of this work and results of work |13j (d(U/w)- 
dependence, the condition of metal-insulator transition, etc); at the same 
time expression fll8|) depends on polar state concentration similarly to re- 
spective expression in work flT3|. 
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3 Energy gap and polar states concentration 



The energy gap (difference of energies between bottom of the upper and top 
of the lower Hubbard bands) is given by 



AE = E d (-w) - E h (w) = -2w(l - 2d) + VU 2 + w 2 , (19) 

(where w = z\t\ is the halfwidth of uncorrelated electron band, z is the 
number of nearest neighbors to a site). Expression (|19"D describes the vanish- 
ing of the energy gap in the spectrum of paramagnetic insulator at critical 
value (2) c when the halfbandwidth w increase (under pressure). Depen- 
dence of AE on temperature can lead to the transition from metallic to 
insulator state with increase of temperature (in this connection note the 
transition at increasing temperatures from the state of paramagnetic metal 
to the paramagnetic insulator state in the systems NiS2- x Sea;, (Vi_ x Cr-^Oa 
and Yx^Ca^TiOa). 

For the calculation of polar states concentration we use function flTT|). At 
T = and rectangular density of states the concentration of polar states is 

, 1 U /l-4cA 

^=7 + ^ ln h^l (20) 



4 8w \ 3 - Ad 



if (2) < m c and 




(21) 



if (2) > (^) c . At T = we have (^) c = 1.672 

u: a- y if 



The dependence d{— ) given by Eqs.(pO|)-(pTD is plotted on Fig. 1. One 
can see that in the point (— ) c the slope of d(^)— dependence changes; the 
concentration of doublons vanishes at — — > oo. Our result for d(— ) in reerion 
of MIT is in good agreement with result of papers |], ^| obtained in the limit 



of infinite dimensions (Fig. 2). The parameter U is normalized by averaged 
band energy in absence of correlation eo- 

In Fig. 3 the dependences of polar states concentration on parameter 
— at different temperatures are presented. Note the important difference 
(see Fig. 4) of the dependence of d on temperature from result of papers || 
E3[: we found that at any temperature polar states concentration increases 
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monotonically with increasing temperature at the fixed value of — when 
respective dependence in ||, 23] has a minimum. 

The dependence of ^ on parameter ^ at zero temperature is plotted 
in Fig. 5. It is important to note that in the point of gap disappearence 
d 7^ in contrast to the previously obtained result [pL3[| . At increasing — 
the energy gap width increases (the negative values of AE correspond to the 
overlapping of the subbands). For comparison on Fig. 5 results of approxi- 
mation "Hubbard-I" J7j is also plotted. In the point of energy gap vanishing 



1.672 what is very close to result of "Hubbard-III" approximation [g2 



At increase of temperature in metallic state the overlapping of subbands 
decreases and temperature induced transition from metallic to insulating 
state can occur at some values of parameter — (Fig. 6). The obtained re- 
sults allows us to draw the (w/U,T) phase diagram of the model (Fig. 7). 
This phase diagram can explain the experimentally observed transitions from 
metallic to insulating state with increase of temperature and from insulat- 
ing to metallic state with increase of bandwidth (under external pressure) in 
paramagnetic state. 



4 Ground state energy 

The ground state energy of the model 



^ = jf(L*H^ia) + Ud 7 (22) 
ija 

calculated using single particle Green function ( [T7|) and expressions (pO|)-(pl|) 
for the concentration of polar states has the form: 

N~ 2 + A {l + M) 2^4(1-2^-1 { ' 



if m < (Z) c and 



^ = -\VWT^ + 2U{\-d) (24) 



if (^) > (— ) c . In Fig. 8 the dependence of the ground state energy on 
parameter ^ given by Eqs.(|23|)-(|2"4"|) is compared with the exact result, found 
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in one-dimensional case |24|. The upper and lower bounds on ground state 



energy in one-dimensional case found in paper |25j are also shown. Our result 
for the ground state energy in metallic state lies slightly lower than exact one 
and in insulator state fits the exact ground state energy very well. 

In Fig. 9 our plot of the ground state energy is compared with the best up- 



per and lower bounds on ground state energy in infinite-dimensional case [26 



In Fig. 10 we have the comparison with bounds on ground state energy for 



three-dimensional simple cubic lattice obtained in paper [25]. In Figs. 8-10 



the ground state energy per electron is normalized by averaged band energy in 
absence of correlation e ; in considered case and rectangular density of states 
Eq = — f . Figs. 8-10 show that our result present a good approximation for 
the ground state energy of the system. In Fig. 11 we plot our result for the 
kinetic part of ground state energy. This plot describes the same behavior of 
kinetic energy of electrons with change of correlation strength in paramag- 
netic state as respective result of work |23[: in metallic state absolute value 
of kinetic energy decrease rapidly due to rapid decrease of doublon (hole) 
concentration. In insulating state absolute value of kinetic energy decrease 
slowly what in the approximation of effective Hamiltonian (obtained for the 
case 7f <C 1) is equivalent to the interaction of local magnetic moments. 



5 Conclusions 

In this paper we have studied electron correlations in narrow energy bands 
using recently proposed approximation fill . We assume that state of the 



narrow-band system is paramagnetic insulator or paramagnetic metal. The 
single-particle Green function and energy spectrum dependent on model pa- 
rameters and on polar states concentration (thus on temperature) have been 
found in paramagnetic state at half-filling (n = 1). The obtained expression 
for energy gap allows to describe MIT at changes of bandwidth (under ex- 
ternal pressure) or temperature. The comparison of calculated ground state 
energy with results of other approximations and the exact result found in 
one-dimensional case shows that the used method is a good approximation 
for the model under consideration. The obtained phase diagram of the model 
can explain the transitions from paramagnetic metal state to paramagnetic 
insulator state at increase of temperature and the paramagnetic insulator 
- paramagnetic metal transitions under external pressure observed in the 



9 



systems NiS 2 - a: Se a; , (Vi- x Ct x ) 2 3 and Yi-^Ca^TiOs. 

It is worthwhile to note that approximation used in this paper can be 
generalized to describe effects of antiferromagnetic ordering. Such a general- 
ization will be considered in subsequent paper. 
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Figure captions 

Fig.l The dependence of doublon concentration d on U/w at zero tem- 
perature. 

Fig. 2 The comparison of d{U/w) dependences: solid line - our result, 
dashed line - iterative-perturbative theory [9,23], circles - QMC method [23]. 

Fig. 3 The dependences of doublon concentration d on U/w at different 
temperatures: lower curve corresponds to kT/w = 0.16, middle curve corre- 
sponds to kT/w = 0.08, upper curve corresponds to kT/w = 0. 

Fig. 4 The dependences of doublon concentration d on temperature at 
different U/w: values of U/w from down to up are 3, 2, 1.5, 1, 0.5, 

Fig. 5 The dependences of energy gap width on U/w: "Hubbard-I" ap- 
proximation (upper curve), our result (middle curve), approximation [13] 
(lower curve). 

Fig. 6 The dependences of energy gap width on temperature at different 
U/w: values of U/w from down to up are 0.5, 1.2, 1.5 

Fig. 7 The obtained (kT,w/U) -phase diagram of the model. 

Fig. 8 The comparison of ground state energies in one-dimensional case: 
dashed curves correspond to upper and lower bounds given by Langer and 
Mattis [25], upper solid curve corresponds to exact ground state (Lieb and 
Wu [24]) lower solid curve corresponds to result of this paper. 

Fig. 9 The ground state energy found in this paper (upper curve), best 
upper (middle curve) and lower (lower curve) bounds on ground state energy 
in infinite-dimensional case. 

Fig. 10 The upper (upper curve) and lower (lower curve) bounds on ground 
state energy in three-dimensional case [25] and the ground state energy found 
in this paper (middle curve). 

Fig. 11 The kinetic part of ground state energy as a function of U/w. 
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